Vortex lines in the three-dimensional XY model with random phase shifts 



Mai Suan Li 

Faculty of Engineering and Design, Kyoto Institute of Technology, Matsugasaki, Sakyo-ku, Kyoto 606 Japan 

Institute of Physics, 02-668 Warsaw, Poland 



or 

< 



Raymond and Beve 



Thomas Nattermann 
Institut fur Theoretische Physik, Universitat zu Koln, 50931 Koln, Germany 

Hciko Rieger 

HLRZ c/o Forschungszentrum Jiilich, 52425 Jiilich, Germany 
Moshe Schwartz 

Sackler Faculty of Exact Sciences, School of Physics and Astronomy, Tel Aviv University, Ramat Aviv, 

Tel Aviv 69978, Israel 
(February 1, 2008) 



The stability of the ordered phase of the three-dimensional XY-model with random phase shifts 
is studied by considering the roughening of a single stretched vortex line due to the disorder. It is 
shown that the vortex line may be described by a directed polymer Hamiltonian with an effective 
random potential that is long range correlated. A Flory argument estimates the roughness exponent 
to £ = 3/4 and the energy fluctuation exponent to u) = 1/2, thus fulfilling the scaling relation 
u) = 2£ — 1. The Schwartz-Edwards method as well as a numerical integration of the corresponding 
Burger's equation confirm this result. Since £ < 1 the ordered phase of the original XY-model is 
stable. 
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PACS numbers: 05.70Jk, 64.60Fr, 64.70Pf 
I. INTRODUCTION 

Systems with a continuous symmetry of the order pa- 
rameter arc particularly susceptible to the influence of 
frozen-in disorderu. In this paper we investigate the sta- 
bility of the ordered phase of the three-dimensional XY- 
model with random phase shifts with respect to topolog- 
ical defects. 

The Hamiltonian of this model may be written in the 
following form 



H = — cos( 

(ij) 



-Oj-Aij) 



(1) 



where J is an effective coupling, the phase variables 
0i G [0, 2ir] are placed on a 3-dimensional hypercubic lat- 
tice and the sum runs over all nearest neighbor pairs (ij) . 
The variables A^ are quenched random phase shifts (or 
random gauges) on the bonds connecting nearest neigh- 
bors. For simplicity we assume, that the Aij on differ- 
ent bonds are uncorrclatcd and gaussian distributed with 
mean zero and variance a. Below it will be useful to go 
over to the continuum description in which A^ is replaced 

by the field A(r). 

Model (Q) describes XY-magnets with random 
Dzyaloshinskii-Moriya interactiono. Other realizations of 
this model are-SD-Josephson-junction arrays with posi- 
tional disordem and vortex glasses-! n 

In the case of the so-called gauge glass modela one 
assumes Aij to be uniformly distributed between and 
2ir but we expect that our model with gaussian disorder 



is equivalent to the gauge glass model when a is large 
enough. 

In d — 2 dimensions the model (0) shows a low tem- 
perature weak-disorder phase with quasi long range or- 
der (QLRO) and, for cr™D a vortex-driven Kosterlitz- 
Thouless-like transition mm to the disordered phase. A 
finite amount of disorder shifts the transition to lower 
temperatures and reduces the universal jump in the spin- 
spin correlation exponent from r/ — 1/4 for a to 
r] = 1/16 at <j c = 7r/8. As it was shown recently QtD, tt 
transition is not reentrant, contrary to earlier findings I 
but in agreement with results of Nishimori for a model 
with slightly different correlations of the random phase 

shifts la. 

In d — 3 dimensions, it is easy to show that spin 
wave excitations, which couple only to V • A, do not de- 
stroy true long range at all temperatures and disorder 
strengths. To address the question, whether this picture 
is qualitatively changed even for weak disorder and low 
temperatures by allowing vortex configurations we con- 
sider a single strechted vortex line in the presence of the 
disorder field A. Such a vortex line can be forced into the 
system e.g. by appropriate boundary conditions. As long 
as the vortex line remains self-affine with a roughness ex- 
ponent £ < 1, which implies (for weak disorder and low 
T) a finite line tension, we conclude, that the ordered 
phase remains stable. It turns out, that this problem 
can be mapped approximatively onto a directed polymer 
problem with long range correlated disorder. We find 
both from analytical and numerical calculations, that the 
roughness exponent is indeed about 0.75, such that the 
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ordered phase remaips stable. This is in agreement with 
results of Nishimori H for a slightly different probability 
weight for the disorder. 



(H) is decomposed into two independent spin wave and 
vortex parts. In what follows we will be interested only 
in the vortex part Tl v which may be written in the form 



II. SPIN WAVES AND VORTICES 

In order to separate between spin wave and vortex de- 
grees of freedom we start from the continuum description 
of model (Q) 

H = J - J d 3 r\V9(f) - A(f)\ 2 , (2) 

The quenched vector field A(r) is assumed to be Gaussian 
distributed with zero average and correlations (a, ft = 1, 
2, 3) 

(A a (r)Ap{f>)) = a5 a0 5(f-P), (3) 

where (• • •) denotes the disorder average. 

The original model (||) is periodic in — Oj with the 
periodicity 2tt. In order to preserve this periodicity in the 
continuum version (^|) we have to allow for singularities 
along which V# jumps by 2ir. These singular surfaces 
are bounded by vortex lines I which are characterized 
by their topological charge m; and their position vector 
Ri(s). It is convenient to decompose the 6*-fleld into a 
spin- wave part V9 SW and a vortex part V9 V . The spin 
wave part is vortex free. The vortex part of the vector 
field V#v is defined by the saddle point equation 

V • (V0„ - 1(f)) = (4) 

and 

V x (V9 V ) = 2-Km(r) , (5) 

where rn(f) denotes the vortex density field which is non- 
zero only along the singular lines Ri(S) 

m(r) = Y / m l j ' ds^^S™ (r - R t (s)) (6) 

Eqs. (||) and @ can be solved easily by introducing a 
vectorpotential a(r ) , Vxa(r) = V9 v —A(f). The solution 
is 

V9 v (f) = J d 3 r'G(r- P){-V'(V' ■ A(r')) 

+ 2ttV' x m(r')} (7) 

with V-m = 0. Here V denotes a derivative with respect 
to the primed variable r' , x represents a vector product. 
The Green function G(f) satisfies V 2 G(r) = — 8(f). In 
three dimensions G(f) takes the form G(f) = (Airr)^ 1 

Representing 9(f) as a sum of a vortex part 9 v (f) and 
a spin- wave part 9 sw (f) one can show that Hamiltonian 
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H vv = £ / / dW^^l , (8) 
^ J J \r — r'\ 

n vd = J - [ ffaV^-^/'^ . 

2 J J \r-r'\ 

The lattice version of such a vortex Hamiltonian for the 
gauge glass has been derived recentlyu. Obviously, TL VV 
and Ti v d correspond to the vortex-vortex and vortex- 
disorder interactions respectively. As in d = 2 dimen- 
sions, the vortices couple only to the V x A(f) — 2irQ(r), 
which can be considered as a quenched random vector 
charge field with V • Q = 0. 

As follows from (3), {Q{r)) — and 

{Q a (r)Qp(r')) = ^-[d a dp - 5 a0 V 2 ]5(r- ?) (9) 

Eqs. (8) and (9) would be the starting point for the 
statistical treatment of model (1). In the partition func- 
tion we had to integrate over all possible configurations of 
vortex loops and vortex lines spanning the system. How- 
ever, this task is much too difficult and remains to be 
done even for the pure system. 

Instead, we follow here a much more modest approach 
and consider the case of very strong dilution of vortex 
lines. In order to test the stability of the vortex free 
state with respect to vortex formation, it is indeed suf- 
ficient to consider an isolated large vortex loop. With- 
out disorder, such a vortex loop costs an energy of the 
order ~ TjqLXyiL (see below) if L is the radius of the 
loop So ~ Jm 2 denotes the bare line tension of vortex 
line. The configurational entropy of the loop is also of 
the order L In L and hence will produce a negative free 
energy only at sufficiently high temperatures. At these 
temperature the system is then disordered, the state of 
the system is characterized by multiply entangled vortex 
lines and loops. 

However, at low temperatures this mechanism does not 
work. A possible source for the condensation of vortex 
loops here is the disorder. To check this possibility, we 
consider in the next section a single stretched vortex line 
which we allow to become rough under the influene of the 
disorder. As long as the typical transverse distortion u ~ 

of a piece of length t of the vortex line is characterized 
by a roughness exponent £ < 1, the contribution ST, t 
of the disorder to the total vortex line tension, S, from 
distortions on this length scale is of the order — Tjqu 2 /t 2 cx 
__^2(c-i) anc j nence small for large t. 

Summation over the contributions <5E t from all length 
scales b n a (a being the microscopic length scale or small 
scale lattice cut-off and b the usual renormalization fac- 
tor) between n min = and n max = ln(L/a)/ln6 yields 
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a finite value for <5£totai which can be made arbitrarily 
small for decreasing disorder strength. Thus the line ten- 
sion £ = So + <5Etotai remains positive and the system is 
stable with respect to the condensation of vortex lines. 

However, for C = 1 the disorder contributions to the 
line tension arc independent of t and hence the energy per 
unit length will vanish on a sufficiently large length scale 
t c . As a result, vortex loops of size L > t c will appear 
spontaneously and hence destroy the ordered phase. In 
this way we traced back the existence of a vortex free low 
temperature phase to the determination of the roughness 
exponent £ of a single stretched vortex line. We expect 
no difference in the case of a vortex loop as long as t c is 
sufficiently large, which can be always acchieved for weak 
disorder. 



III. SINGLE VORTEX LINE HAMILTONION 

We consider in the following only a single stretched 
vortex line with no overhangs, which means that we 
may set R{s) — > R(t) = (ri(t), r 2 (i), t), where r(i) = 
( r i(t), r 2{t j) describes the transverse distortion of the 
vortex line from a straight configuration. Then rfi(r) may 
be parametrized as followsEj 



m(r) 

dR 
~dt 



ln / dt^5^(r- R(t)) . 
at 



dr\ <9r 2 



(10) 



H vv decribes now the elastic self-interaction of the vortex 
line H vv — ► H e i- 



H e 



dR(t)dR(t') 



Jtt 



dr dr 
di'di* 



dR{t)dR(t') 
\R(t)-R{t')\ 

1 ^1 dtdt' . (11) 



Expanding H m in Eq.(||) up to quadratic terms in 
r and dr/dt and omitting an irrelevant constant corre- 
sponding to the energy of straight lines one has 



He 



Jtt 



dtdt' ,dr dr 



\t-f\ dt dt> 
(rft)-r(Q) 2 
2\t~t'\ 2 



0(r 4 



(12) 



It is convenient to go over to the Fourier transform r(i) 



Jo l ' 



H e l — 



Jm 1 



du> 



w 2 /(aw)r w .r_ u 



2 J (2tt) 

where a denotes a small scale lattice cut-off and 



(13) 



dx 



[(1 



1 



1 



lne-- (14) 
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for e -> 0. 

Next we consider the correlations of the vortex disorder 
interaction. For this purpose we rewrite H v d as 



H vd = dR- V(R) 



where 



V(R) = Jirm / d A r 



,3„ Qif) 



\R-f\ 



(15) 



(16) 



Clearly (V) = and 

, - s , TxJ 2 m 2 1 

(V a (R)Vp(R>)) = a - ift ~ {S a0 + 



2 \R-R', 
(R a — R' a )(Rf3 — R'p) 

(R - R') 2 



} (17) 



Since the correlator (|l7| ) always appears under the inte- 
gral over R, it is convenient to rewrite the right hand side 
©as 



r {Jm) { dRiM f 



\R-R' 



\R-R' 



;) 



The first contribution in this expression leads to terms 
in the correlator of H v d which vanish by assuming pe- 
riodic boundary conditions. Thus we omit it in the fol- 
lowing. For small gradients |^| <C 1, H v d can finally be 
rewritten as 



H v d = / dtV(r(t),t) 



(18) 



with 



(V(r,t)V(r',f)> = a{Jmf\{v - r') 2 + (t - t') 2 ^ 2 

(19) 

Thus, the random potential V(R) which interacts with 
the vortex line is long range correlated. 

IV. FLORY ARGUMENTS AND BEYOND 

In the following we want to estimate the conditions un- 
der which the random potential can destabilize the vor- 
tex line. We begin with a straight line of length t, which 
has an elastic energy E e i — ^-m 2 t\nt/ a as follows easily 
from ([ll]) . The typical fluctuations of the disorder energy 
follows from ( |l7| ) as E d i S = ±mJy / 7ro : (t lnt/a) 1 / 2 . Since 
rare fluctuations increase Edi S only by a factor (lni)' 1 / 2 - 1 , 
the system is always stable with respect to the formation 
of straight vortex lines. 

Next we allow for displacement r(t) from the straight 
configuration such that |^|| <C 1. If on scale t the vortex 
line is self-affine with u — (r 2 ) 1 / 2 ~ t^ with C, < 1, we can 
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use a Flory-argument to estimate C- The elastic energy 
is then increased by AE e i w J " 7r ^~ 1 u 2 lnt/a which has 
to be compared with ii^is- This yields 



u ~ [(rt 3 (lni/a) 



-111/4 



(20) 



i.e. C — 3/4- Although the Flory argument is rather 
crude, it seems to be safe enough to conclude ( < 1, i.e. 
for weak disorder the vortex line is self-affine and hence 
stable (only £ — 1 would signal an instability with respect 
to vortex generation). 

To go beyond the Flory argument one usually maps 
these types of directed polymer problems to Burger's 
equation with (correlated) noise. However, in the present 
case this is not possible in the strict sense because of 
the long range elastic self-interaction of the vortex line, 
which leads to the factor f(au>) w In -i- in (|l4|). Since 
this corresponds to a logarithmic increase of the stiffness 
of the vortex line with the length scale, we neglect the 
cj-dependence of this factor in the following completely 
and set equal to /n. If then the roughness exponent £ is 
less than one, it will be less than one by including f(auj). 
In fact, it is even safe to assume that the value of C is 
unchanged by replacing f(au>) by /n since for small auj 
f(atu) diverges only logaritmically. 

With this simplification the free energy h(r, t) of a vor- 
tex line of length t and r = r(i) — r(0) obeys the Burger's 
equation with noise 

?- t h(r,t) = vV 2 h+~(Vh) 2 + V(r,t) (21) 

where v = T/2Jm 2 fo and A = 1/Jm 2 f . 

Besides of the roughness exponent £ there is a second 
exponent w, which describes the sample to sample fluc- 
tuations (h 2 (r,t)) — (Mr,t)) 2 ~ t 2u of the free energy. 
Since the correlator (|19| ) is non-local in time, the validity 
of the scaling relation oj = 2Q — 1 is not guaranteed!!^ . 

If we give h(r, t) another interpretation, namely that of 
the height profile of a growing interface in a cq-wioving 
frame, eq. ( pl| ) is known as the KPZ-equationli3. The 
roughness exponent a, the growth exponent (3 and the 
dynamical exponent z of this surface are related to the 
directed polymer exponents by 



C = (3/a = l/z 
lu = (3 = a/z . 



(22) 



We have first attempted to use the standard one-loop 
renownalization group treatment for the Burger's equa- 
tioiifcy. Unfortunately, this approach does not yield (as 
happens also for the case of uncorrelated noise in three 
dimensions) stable fixed points from which one could cal- 
culate critical exponents. We study the Burger's equa- 
tion with the correlated noise (eq.21) therefore here by 
a method that proved useful for the same problem with 
uncorrelated noise L3. 

Our first step here is to describe our system as a system 
of dynamical variables affected by noise that is uncorre- 
lated in time. This is necessary in order to go over from 



the Langevin-like Burger's equation to the Fokker Planck 
equation, needed in the method mentioned above. 

This step is achieved by considering the V to be also 
dynamical variables coupled to a noise £ that is uncorre- 
lated in space and time. 

In Fourier components we write eq. ( ^l"l) in the form 



dhq 

dt 



-vtfK - — 2 1 ■ (q - 1)^., - Vq(t) 



1 



(23) 



where f2 is the transverse area of the system(to be taken 
eventually to infinity) and add to it an equation for V q (t) 



BV 

-gf = -|q|V q + U*) 



(24) 



where (f q (t)) = and (£ q (t)£- q (t')> = 2D5(t - t') and 
2D is short for 47t<t(Jto 2 ). It is easily verified that Vq(t) 
has the requiered correlations given by eq.(19). 

We write now the Fokker-Planck equation for the joint 
probability density to have a given h and V configuration, 

+ ^E 1 -(q- 1 )^-' + ^)} F ( 25 ) 

We will be interested in the steady state averages: 
{hqh-q) s , {hqV-q} s , and (VqV-q) s and in some char- 
acteristic frequencies to be defined later. We assume 
now that the exact values of (hqh—q) s ,{hqV—q) s and 
(VqV-q) s are known and given by X" 1 ^" 1 and A" 1 . 
In order to calculate such quantities, we need some form 
of a perturbation expansion in which the Fokker-Planck 
operator O, acting on P on the right hand side of eq.(24), 
will be broken into two parts: one part Oq that is simple 
enough and another part O — Oq, that is small enough. 
We may expect that if we choose Oq in such a way, that it 
already gives the exact result for the three steady state 
averages, the corrections to those quantities in pertur- 
bation theory, are going to be small giving sense to the 
whole expansion. 

Our choise for Oq is 



On 



,(1)J_ [A -1 9 

q «t/ 1 q 



q 



dV-q 



V a 



(2) 







d 



q dhq 1 q dk-q 



,(1) 



,(2) 



][A. 



q dh~qdV7q\ 



(26) 



It is easily verified, that regardless of the choice of the 



u)'s, Oo produces {hqh- 
and (VqU_ q ) s = A q !. 



q/s 



X~\(hqV-q) s = A" 1 
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A direct inspection of eq. (24), shows that the form 
of the V dependent part of O is identical to the form of 
the correspoding part of Oq. Therefore we must identify 
Aq 1 ^ and 4 X) = |q|. 

Next we calculate in perturbation theory (/i q /i- q ) s and 
{hqV-q) s . We find 



(hqh-q) s 

(hqV-q) s 



(2) n 



A-' + CziX- 1 ^- 1 ,^} 



(27) 



Since A q 1 and A q 1 are assumed to be exact, we get 

(2) 

now, for a given set {w q }, two equations C\ — C2 = 
expressing the fact that the perturbation expansion does 
not change the zero order value. Actually it could be 

expected that for any choice of positive ojq s, the above 
equations for the A q s and A q s(Ci = Ci = 0), would 
give the exact result. It is clear however, that even for 

(2)' 

statics, we may use the choice of the ojq s in order to 
control the smallness of the additional terms and thus we 
will be able to obtain our X^s and A q s from low order 
perturbation. Assuming that such a choice is possible 
we write C\ and Ci to second order to obtain the two 
equations: 



0= - 



•4 2) ]a- 



r 2 (2)i2 

WQ ] A-l 

- 7T\ 



}D 



|q|[<'+< 



D 
|q| 



(28) 



and 
= - [2vq 2 

+ N 2 -4 2) ][ -fa + -mb^r ]A- X 



2 - { ^\Xq 1 



+ 



|qlK 1J +<- 



2q • (1 - c^X-^X- 1 



(2) 

We fix now ojq to be the characteristic frequency as- 
bions 

(fe q fe_ q ) s 



(29) 



sociated with the correlations of hq at different times 



,(2) 



q f °°(h.q(0)hq(t))dt 
where (A(O)B(i)) is defined for t > to be 



(30) 



(A(0)i?(t)) - / A{h' v V£P a {h' v V£P{h' v hq, y q ; t} 
Bih^VJDh'^DV^DhqDVq . 

(31) 

P s is the static (steady state) distribution and 
P{h' , F q ; hq, Vq, t} is the solution of the Foker-Planck 



equation for initial conditions of absolute certainty that 
the value of hq is /i q and V q is (II q #(/i q — h'q)6(Vq— Vq) 
is the initial distribution). Our lowest order Fokker- 
Planck operator, Oo, is thus chosen also to give the cor- 
rect characteristic frequency. We calculate now the char- 
acteristic frequency to second order and demand again 
that the correction to the zero order term is zero. This 
yields 

= 



2(2tt) 



+ 2(2tt)2 / d 2 l- 



'l'(q-l) 



[2q • (1 - q)Ay r A q 1 + I- (q - l)X{~ L X~^ l ] 



Apart from the last term, eq. (| 
corresponding equation ir£3. 
We assume now that for small q 



(32) 

3) is identical to the 



(2) 
^q 

A„ 



cx q z + 
oc q Tl + 



Aq oc q v 



corrections 
corrections 
corrections 



(33) 



The exponents z and Ti are related to the exponents £ 
and ui by the relations C — X/z and u> = (Ti — 2)/2z. We 
solve equations (^8|), (^9|), and (^2|) to leading order and 
obtain the folowing results. The scaling relation 



(34) 



is obeyed . (The last term in eq. ( p2| ) is negligible compared 
to the dominant terms for small q). This is the familiar 
scaling relation 



uj = 2C - 1 



(35) 



Eqs.(p8|) and (29) allow now for a simple power counting 
solution (In contrast to ref.13 where a power counting 
solution was impossible). 
We obtain 



r,-» * = 1, r a = 2 



(36) 



It is interesting to note that the result for I^ is not 
affected by A being not equal to zero and this is because 
the last term in eq.(p8[) is negligible, for small q's, como- 
pared to the dominant terms. Consideration of higher 
order terms in the expansion, (in particular, the most 
dominant q dependence in each order) seems to suggest 
that this results will still hold to any order of pertur- 
bation theory. In terms of the more familiar directed 
polymer exponents we find 



t 3 1 
C = — uj = — 

s 4 2 
so that the Flory result still holds. 



(37) 
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V. NUMERICAL RESULTS 

In this section we will obtain the roughness exponent 
for our model by integrating the (2+l)-dimensional KPZ 
equation ( pl| ) numerically. The simplest quantity to in- 
vestigate is the surface width 



W = (h 2 - h 



■2' 1/2 



(38) 



where the bar and angular brackets (. . .) denote a spatial 
and noise averages respectively. One expects for Wl the 
following scaling formEJ 



W = L a f(t/L z ) 



(39) 



where f(x) ~ x@ with (3 = a/z for x <C 1 and f(x) ~ 
const, for x 3> 1; /3, a and z are the growth, roughness 
and dynamic exponent of the interface respectively, c.f. 

In analogy to (^) the height-height correlation func- 
tion 



C L (r, t) = ([h(r + ?, t + t')- h(P, t')] 2 ) 1 ' 2 
should scale as 

C L (r,t) = ^gir/f/^L/t 1 /*) . 
We do not take the scaling relation 
a + z = 2 , 



(40) 
(41) 
(42) 



which is equivalent to (|35|), for granted in our investiga- 
tion. Thus, to obtain a full set of the critical exponents 
we have to calculate at least two oL-theta.. 

Different theoretical approachesOE3ii3 lead to differ- 
ent dependences of the critical exponents on p in (1+1)- 
dimensions. To check theoretical predictions the numer- 
ical .simulations have been carried out for the KPZ equa- 
tion 11, the ballistic depositionLjCDtp, the directed poly- 
merEZI and the solid-on-solid modeled. There is still a con- 
troversy between the simulation results. For example, the 
results of Amar et alxB obtained for the ballistic deposi- 
tion and restricted solid-oa,-solid models agree with the 
prediction of Medina et alM but conflict with the predic- 
tion of ZhangtJ. Numerical studieso of the effect of the 
long-range spatially correlated noise on the KPZ equa- 
tion and the related directed polymer problem, on the 
other hand, give a good-agreement with the prediction of 
Hentschel and FamilyEJ. Temporally correlated noise in 
the absence of spatial poise correlations has also been in- 
vestigated numerically 2 ^ for the (l+l)-dimensional case. 

Theoretical results for (2+l)-dimensional models f*4thi 
correlated noise are still lacking. Meakin and JullienliSEil 
introduced a hopping model of ballistic deposition, in 
which particles were deposited on the growing surface fol- 
lowing a Levy flight distribution such that the distance x 
(along the surface) from the previous site was calculated 
as x = r^ 1 ^ where r is a random number between zero 



and 1_ Equating the exponent / to 2 as in ( |l9| ) led to re- 
sultsa for the exponents a and (3 which were roughly in 
agreement with the prediction of Medina et al. for (1+1) 
dimensions. A weak dependence of these exponents on / 
has been found in (2+l)-dimensional modelcil. 

In the simulations of Meakin and Jullien the link be- 
tween the deposition process and the noise was not ob- 
vious. Thus the scaling behavior in (2+l)-dimensional 
interface growth and the directed polymer model with 
spatially correlated noise is at the present time unclear. 
We consider this problem solving the KPZ equation nu- 
merically. The spatial derivatives in ( pl| ) are discretized 
using standard forward-backward differences on a hyper- 
cubic grid with a lattice constant Ax. Eq. (|2l]) is in- 
tegrated by the Euler algorithm with time increments 
At. Denoting the grid points by n and the basic vectors 
characterizing the siHjface by ei, . . . , e<j we arrive at the 
discretized equationES 



h(n,t + St) = h(n,i) 

>t ir^d , 

+ hh(n + e i7 t) - h(n- e 4 ,t)] 2 } + V3Atr)(n,i) 



+ Eti { [Hn + St, i) - 2h(n, t) + h(n - e h t)\ 



(43) 



Here one uses dimensionless quantities h — h/ho,x = 
x/xq and t = t/to where the natural units are given by 



v v 
ho = —, t = , xq 



a 2 = 2D/Ax d 



(44) 

Similar to V(r,t) in the Fourier space the renormalized 
noise r](n, t) has the correlation 



(n(k,u})r](k',uj')) 



S(k + k')5(u + u)') 
k 2 + uj 2 



(45) 



Tcj-Gj-eate the spatially correlated noise we follow Peng et 
alS-H. We first generate a standard white (or Gaussian 
uncorrelated) noise T]a{n,t), then carry out the Fourier 
transformation for spatial and temporal variables to ob- 
tain Tfo(k,u>). We define 



V (q,co) = (fc 2 +^ 2 )- 1 / 2 % (fc,^) 



(46) 



The noise ry(n, t) is obtained by Fourier transforming 
i](k,uj) back into the real space. It is easy to check that 
rj{k,ui) obained by this way satisfies ([45|). 

In our simulations we chose Ax, v and a to be the same 
as in Ref.E3, namely Ax = 1, v = 0.5 and a — 0.1. To be 
sure that we are in a strong-coupling regime we chose A = 
v600. The data we show in the following are averaged 
over 10000 samples (i.e. different noise realizations) for 
the smallest system size (L = 16) to 50 samples for the 
largest system size (L — 256). 

First we determine the exponent P and the dynami- 
cal exponent z from simulations of large system sizes on 
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FIG. 1. Log-log plot of the width Wh(t) vs t for various 
system sizes L. For smaller times the data are roughly in- 
dependent of system size with a slight tendency to a smaller 
slope for increasing L. For larger times (i.e. when t ~ L z ) sig- 
nificant finite size effects appear (i.e. the scaling function w in 
( [ill ) is not constant any more. Therefore we put more weight 
on the data points at smaller times in order to determine 0, 
which yields = 0.56 ± 0.02 shown as dotted line. 



short time scales, i.e. we focus on a regime where the 
correlation length £ is still small compared to the system 
size L: 



(47) 



In this regime finite size effects should be negligible, 
which we had to check explicitely by analyzing differ- 
ent system sizes since we do not know z a priori. For the 
time dependence of the roughness we expect the scaling 
form to hold, which can be cast into the form 



W L (t) = t> 3 w(t/L z ), 



(48) 



with w(x) = const, for x — > and w(x) ~ x~ 1 for x — > 
oc. For short times t <C L z one therfore expects Wz,(t) oc 
t&. Fig. 1 shows the time dependence of the width Wl (t) 
for various system size L at short time scales. We find 



= 0.56 ±0.02 



(49) 



This value is much larger than ~ 0.240 for the three- 
dimensional KPZ equation with the white noiseE3, as we 
expect it to be since correlated noise always increases 0. 

In order to obtain the dynamical exponent z from the 
short time behavior we note that for t -C L z the height- 
height correlation function C(r, t) (^) should have the 
scaling form 



C(r,t) = r^g{r/t l l z ), 



(50) 



where we denoted lim( tL -z)^ ^(r/t 1 ^, L/t 1 '*) = 
gir/t 1 /*) and used a = 0z. 

If we take j3 = 0.56 as determined before we are left 
with one fitting parameter by which we should achieve a 




0.001 

FIG. 2. Scaling plot for the height- height correlation C(r J t) 
function. We fixed /3 to be 0.56, as determined before (149) 
and varied z in order to achieve the best data collaps, which 
yiels z = 1.25 ± 0.05. 

data collapse when plotting r~P z C(r, t) versus r/t 1 ^. In 
fig. 2 we show the result of this procedure, which yields 



z = 1.25 ±0.05 



(51) 



Si 

-J 




FIG. 3. Scaling plot of Wi(t) according to the scaling form 
©. The choice a = 0.71 ± 0.01 and z = 1.25 ± 0.03 yields 
the best data collaps, which is shown. 

Up to now we have restricted ourself to a regime, where 
finite size effects are negligible. In order to get an in- 
dependent estimate for the exponents reported above 
we performed also a simulation for much longer times 
and small enough system size, so that the correlation 
length becomes indeed comparable to the system size (i.e. 
t ~ L z ). Using the finite size scaling form (|39| ) we deter- 
mine the exponents a and z. In fig. 3 we show a scaling 
plot of Wl (t) which yields 



a = 0.71 ±0.01 



(52) 



and z = 1.25±0.03, from which one concludes — a/z = 
0.56 ± 0.03. These values agree very well with the values 
reported in (49) and ( j5l| ) for the short time simulations. 
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Using scaling relation ( p^ ) we have 



u) = 0.56 ±0.02 and C = 0.78 ± 0.04 (53) 



for the vortex lines in our model. The value of £ is larger 
than the corresponding value ( w 0.62 for the directed 
polymer paths with uncorrelated noise in three dimen- 
sions. Our numerical result is compatible within the er- 
ror bars with the crude estimate £ = 3/4 obtained by the 
Flory argument and by the more elaborate calculation 
in the last section. The numerical value of u is slightly 
larger than the calculated u> — 1/2. 

Note that within the error bars relation fli^ ) is still 
valid even in the presence of the temporal correlation of 
noise. This relation may hold true only for some subset 
of all possible correlators (45). At present we have to 
leave this question open. 



VI. SUMMARY AND CONCLUSION 

In the present paper we have considered the stabil- 
ity of the ordered phase of the XY-model with random 
phase shifts. After decomposing the Hamiltonian into a 
spin- wave and a vortex part, we have considered in par- 
ticular the roughening of a single stretched vortex line 
due to the disorder. It turned out that the effective ran- 
dom potential acting on the vortex line is long range 
correlated. Using a Flory argument and the Schwartz- 
Edwards-method, we have determined the roughness ex- 
ponent Q = 3/4 and the energy fluctuation exponent 
u> = 1/2, which fulfill the scaling relation io — 2£ — 1. 
These findings have been confirmed by integrating nu- 
merically the Burger's equation. Since £ < 1, there will 
be no spontaneous condensation of vortices for weak dis- 
order and hence we conclude, that the ordered phase is 
stable. 

A further interesting application .pf.our result is that 
on the XY-model in a random fieldcjo. From the results 
of Villain and FernandezO for this model without vortices 
one finds a(t) ps a*t d ~ 2 , where a* — O(gEf)- Using this 
t-dependence of a in the Flory-argument one concludes 
( = 1 for 2 < d < 4, apart from a logarithmic correction. 
From this one should expect that one is in a-4narginal 
situation, which deserves further investigationEH 
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